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ABSTRACT 

Context. Stars on the asymptotic giant branch (AGB) lose substantial amounts of matter, to the extent that they are important for the 
chemical evolution of, and dust production in, the universe. The mass loss is believed to increase gradually with age on the AGB, but 
it may also occur in the form of bursts, possibly related to the thermal pulsing phenomenon. Detached, geometrically thin, CO shells 
around carbon stars are good signposts of brief and intense mass ejection. 

Aims. We aim to put further constraints on the physical properties of detached CO shells around AGB stars. 

Methods. The photodissociation of CO and other carbon-bearing species in the shells leads to the possibility of detecting lines from 
neutral carbon. We have therefore searched for the CI(^Pi - ^Pq) line at 492 GHz towards two carbon stars, S Set and R Scl, with 
detached CO shells of different ages, ^ 8000 and 2300 years, respectively. 

Results. The CI(^Pi - ^Pq) line was detected towards R Scl. The line intensity is dominated by emission from the detached shell. The 
detection is at a level consistent with the neutral carbon coming from the full photodissociation of all species except CO, and with 
only limited photoionisation of carbon. The best fit to the observed ^^CO and ^^CO line intensities, assuming a homogeneous shell, 
is obtained for a shell mass of ^ 0.002 M©, a temperature of ^ 100 K, and a CO abundance with respect to H 2 of 10"^. The estimated 
CI/CO abundance ratio is ^0.3 for the best-fit model. However, a number of arguments point in the direction of a clumpy medium, 
and a viable interpretation of the data within such a context is provided. 
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1. Introduction 


It is well established that low- and intermediate-mass stars go 
through extensive mass loss during the evolution on the asymp¬ 
totic giant branch (AGB). Most of this mass loss probably occurs 
over an extended time scale through a wind that appears to in¬ 
crease in strength as the star evolves ( |Habing||l996 ). However, 
there is also ample evidence of highly episodic mass loss. Most 
notable are the geometrically thin, detached CO and dust shells 
around carbon-rich AGB stars. Presently t here are published data 
on seven CO shells ([Schoier et al.||2005|) and 1 7 dust shells of 
this ty pe ( Ramstedt et al. 201 ![ et al.|20f^|Maercker et al. ' 


2014 [Mecina et al.|2014 and references therein) (not ah of the 

dust shells are clearly of the geometrically thin type). We note 
here that some stars with prominent dust shells, e.g. AX Cyg 
and RT Cap ( [Kerschbaum et al.||2010[ |Mecina et al.||2014| , do 
not have detectable CO shells, most likely an effect of CO pho¬ 
todissociation if the shell is old enough. Even though most of the 
mass is not lost during these episodes, they provide an interest¬ 
ing phenomenon which may shed light on the mass-loss mech¬ 
anism of AGB stars, e.g. the mass ejection may be connected to 
a physical process with well-defined variations of stellar charac¬ 
teristics (see below), the line of sight confusion is limited due 
to the thinness of the shells and this allows studies of the small- 
scale structure of the circumstellar medium, and the evolution 


^ Based on observations with the Atacama Pathfinder Experiment 
(APEX) telescope. APEX is a collaboration between the Max-Planck- 
Institut fur Radioastronomie, the European Southern Observatory, and 
the Onsala Space Observatory. 


of the gaseous and dusty components can be compared in some 
detail. 

T he reason for the episodic mas s loss is not definitely iden¬ 
tified. [Olofsson et al.| ( |199Q[ |1993| ) provided good arguments, 
based on the detection statistics of an essentially complete sam¬ 
ple of carbon stars, that the strong modulation of the mass-loss 
rate is driven by the changes that a star goes through during a 
He-shell flash, or thermal pulse. [Schoier et al.| ( |2QQ5| ) provided 
evidence that the interaction of a previous, slower stellar wind 
has an effect on the detached shell characteristics. The mass- 
loss modulation is estimated to be, at least, one order of mag¬ 
nitude, being less than 10“^Moyr“^ before and after the shell 
ejection and larger than 10“^ Mq yr“^ during shell ejection. The 
shell ejection time scale is a few hundred years ( jOlofsson et al. 
|2000[ [Maercker et al.|2012| ). Curiously, no M-stars, which also 
go through thermal pulsing, have been found with detached CO 
and/or dust shells of this type despite sever al searches ([Nyman 
et al.[[ 19921 [ Kerschbaum & Qlofsson[[T999l [Kerschbaum et al. 


et al.[[1992l[Ki 


The physical properties of the detached shells have been 
probed so far using CO radio line emission and through dust in 
either scattered stellar light or thermal emission. The CO lines 
have the advantage of providing information on the kinematics, 
and they can be studied at high angular resolution. Data of the 
latter type show remarkably spherical shells with width/radius 
ratios of < 0.1 ([Lindqvist et al. [1 996 [ [Olofsson et al.[2000tpV[aer- 


cker et al.|2012[ ). The images in scattered light are insensitive to 

any excitation processes ( [Gonzalez Delgado et al.|20OT| [Olofs¬ 
son et al.[2010[ ), and in polarisation mode they provide important 
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geometrical information ([Gonzalez Delgado et al.||2003t |Maer- 


|cker et al.|2010[[20T^|Ramstedt et al.|2Ql 1| ). The images in ther¬ 

mal emission are slightly more sensitive to excitation (the grain 
temperature) and there is more line of sight confusion, but they 
provide good estimates of the dust mass if observed at longer 
wavelengths. It turns out that a determination of the physical 
characteristics (density, temperature, dumpiness, etc.) of a de¬ 
tached shell is difficult using the existing data. 

Hence, there is a need for additional probes, but in the case 
of molecular lines they are hard to find owing to efficient pho¬ 
todissociation. Eventually, CO is also dissociated into C and O, 
though it is much more stable than other species owing to line 
dissociation and self-shielding ( [Marnon et al.|1988] ). Therefore, 
the CI(^Pi - ^Pq) line at 492 GHz is a possible alternative probe 
of the shell characteristics. It is, however, well known that this 
line is not easily detected in AGB and post-AGB objects. There 
are published detections in only one carbon star (CW Le o,|Keene| 
et al.| ( |1993| )), one supergiant (a Ori, [Huggins et al.| ( |1994| )), 


and seven carbon-rich post-AGB objects [CRL618, CRL2688, 
HD56126, HD44179 (Red Rectangle), NGC6720 (Ring nebula), 
NGC7027, NGC7293 (Helix nebula), [Bachiller et al.j (19941); 
|Young| ( [TW7l ); [Young et aLKIW/] ); [Knapp et aLlpOPOl )]. 

In this paper we report on a search for the CI(^Pi - ^Pq) 
line in two detached CO shell sources, R Scl and S Set, where 
the shell ages differ by about a factor of four from the former 
2300 yr) to the latter 8000 yr) ( [Schoier et al.|2005 ). In the 
case of R Scl the search was successful, and we discuss the im¬ 
plications of the Cl detection for our understanding of detached 
CO shells around AGB stars. 


2, Observations 


We have used the APEX 12 m telescop e located at an altitu de 
of about 5100m in the Chilean Andes ( [Glisten et al.[[2006| to 
observe the CI(^Pi - ^Pq) line at 492 GHz. The Swedi sh hetero¬ 
dyne facility instrument APEX-3 ( [Vassilev et al.[200^ , a double 
sideband receiver, was used together with the facility EET spec¬ 
trometer covering about 4 GHz. The observations were made in 
April and June 2014 under programme 0-093.E-9306. 


The observed sources are presented in Table with variable 
type (V), period (P), stellar luminosity (L^), and distance (D) 
for each star. The periods, luminosities, and distances, which are 
estimated by fitting the spectral energy distributions and calcu¬ 
lating the luminosities from the period-luminosity relations, are 
taken from [Gonzalez Delgado et al.[ ([2001[) and [Schoier et al. 
( [2005] ). 


The observations were made in a dual beamswitch mode 
where the source is alternately placed in the signal and the ref¬ 
erence beam, using a beam throw of about 2'. Regular pointing 
checks were made on strong CO line emitters (e.g. in the case 
of R Scl on the object itself) and continuum sources. Typically, 
the pointing was found to be consistent with the pointing model 
within ^ 3". 


The raw spectra are stored in T*-scale and converted to 
main-beam brightness temperature using Tnib = The an¬ 

tenna temperature, T*, is corrected for atmospheric attenuation 
using a calibration unit, and T/mb is the main-beam efficiency. The 
adopted beam efficiency and the EWHM of the main beam at the 
observing frequency of 492 GHz are 0.6 and 13", respectively. 
The uncertainty in the absolute intensity scale is estimated to be 
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Table 1. Source information: variable type (V), period (P), luminosity 
(L^), and distance (D). 


Source 

V 

P 

[days] 

[Lo] 

D 

[pc] 

RScl 

SRb 

370 

4900 

370 

S Set 

SRb 

148 

4300 

400 


Table 2. CI(^Pi - ^Pq) observational results 


Source 

^sys 


^ ^peak 

/ 


[km s“^] 

[km s“^] 

[K] 

[Kkms-i] 

RScl 

S Set 

-19.0 

14.5 

-1.1 0.09 

1.5 

<0.2 


about +20%. The data were reduced by removing a first-order 
polynomial baseline using XSQ 

We note here that frequency shifts of the two isotope 
^^CI(^Pi - ^Po) hyperfine lines with respect to the frequency of 
the ^^CI(^Pi - ^Po) line are small, 0.5 MHz and 3.6 MHz for the 
F = 1/2-1/2 and 3/2-1/2 components, respectively. This means 
that in the circumstellar case where the lines are broad, emission 
from the lines will contribute to the fiux measured for the 
line. However, since the circumstellar ^^CO/^^CO is about 
45 and 20 for S Set and R Scl, respectively, we expect this con¬ 
tribution to be small ([Bergman et al.[[ 19931 [Vlemmings et al. 
[MT3l ). 

3. Observational results 

The observed CI(^Pi - ^Pq) line towards R Scl, binned to a ve¬ 
locity resolution of 1.5kms“^ is presented in Pig. together 
with the CO(/ = 4-3) line at 461 GHz, which is observed with 
roughly the same beam width. The line profile has been fitted to 
the following function. 


T(v) = T(V,y,) 



( 1 ) 


where L'sys is the central velocity, and Vq the gas expansion veloc¬ 
ity of the Cl line emission. A negative a means a double-peaked 
profile. The line intensities are integrated over the line profile. In 
the case of S Set, a 3cr upper limit to the line intensity is esti¬ 
mated. A summary of the results are given in Table The sys¬ 
temic and expansion velocities we obtain for Cl agree well with 
those obtained from the CO lines in the case of R Scl ( [Qlofs-[ 
son et al.][1996[ ). Both lines are double-peaked, but the Cl line 


has a markedly higher peak/centre ratio than the CO(4-3) line 
= compared to -0.4), which indicates that this line is sig¬ 
nificantly more dominated by emission from the detached shell 
(see Sect. HI- 

In addition to these data, we have assembled CO line 
data from the literature to constrain our modelling of the de¬ 
tached shell characteristics. These constitute the ^^CO SEST 
/^(1 -Q) line ([Olofsson et al.[[T996l ), IRAM (1-0) and (2-1) 
lines ([Danilovich et al.|2015[), APEX (2-1) and (3-2) lines (|De 
Beck et al.[[2010[ ), JCMT (3-2) and (4-3) lin es ([Schoier et al. 


2005) and Herschel/HIPI (5-4) and (9-8) lines ( [Danilovich et al. 


2015) towards R Scl, and the SEST / = (1-0), (2-1) and (3-2) 


^ XS is a package developed by R Bergman to reduce 
and analyse single-dish spectra. It is publicly available from 
ftp://yggdrasil.oso.chalmers.se 
























































































H. Olofsson et al.: Cl line emission towards R Scl 




Fig. 1. Observed CI(^Pi - ^Pq) (top) and CO(/ = 4-3) (bottom) lines 
towards R Scl. 

Table 3. Detached shell: gas parameters 


Source 

t^sys 

Vq 

Vi 

Rs 

ARs 


[kms“^] 

[km s“^] 

[km s"^] 

[cm] 

[cm] 

RScl 

-18.9 

14.3 

1.0 

1x10'^ 

7x10^5 

SSet 

15.0 

16.5 

0.8 

4x10'^ 

1x10*® 


lines ( [Olofsson et al.|1996| ) and APEX (3-2) line ( jPe Beck et al.j 
|2010| ) towards S Set. In addition, as further constraints on the 
R Scl modelling (primarily on the opt ical depths) we have in - 


cluded the IRAM / = (1-0) line ([Danilovich et al.||2015|, 
SEST (1- 0) and (2-1) lines ([Qlofs son et al.||1996D and APEX 
(3-2) line ( jVlemmings et al.|2013 ). 


4. Radiative transfer 


We have used a radiative transfer code based on the accelerated 
lambda iteration method (see e.g. jMaercker et al.|2008| ) to derive 
quantitative results. As the geometrical (radius, and width, 
ARs) and kinematical (systemic velocity, Psys, and expansion ve¬ 
locity, Vq) p roperties of the detach ed sh ells we adopt the result s 
obtained by jSchoier et al.j ( |2005|) and [Maercker et al.j p012| ). 
These are summarised in Table l3J The detached shells are as¬ 
sumed to be homogenous in density and kinetic temperature. The 
turbulent velocity width, is chosen to be 1 and 0.8 kms“^ for 
R Scl and S Set, respectively. 

In the and excitation analysis, we included ra¬ 
diative transitions between rotational energy levels up to /= 14 
for the vibrational ground state (radiative excitation by a cen¬ 
tral source is negligible in the detached shells). The energy lev¬ 
els, transition frequencies, and Einstein A coefficients were taken 
from the Cologne database for molecular spectroscopy [CDMS, 
Muller et al.j ( 2005| )] . Collisional rates were taken from [Yangj 
et al. ( |2010| ), who calculated them separately for CO collisions 
with para- and ortho-H 2 . The rates were weighed together as¬ 
suming an ortho-/para-H 2 ratio of 3. They cover temperatures 
from 2 to 3000 K. 


Table 4. Integrated CO line intensities for R Scl in the velocity ranges 
11.1 < I P-Psys I ^ 15.3 km s"^ 


Line 

Tel. 

Beam 

["] 

/obs" 

[Kkms-*] 

ref.® 

‘'*co(y=i-o) 

IRAM 

21 

21.6 + 0.3 

1) 

*2co(y=i-o) 

SEST 

45 

4.6 ±0.1 

2) 

‘2C0(7 = 2-1) 

APEX 

27 

28.2 + 0.2 

3) 

‘2C0(7 = 2-1) 

IRAM 

11 

14.2 + 0.1 

1) 

‘2CO(y = 3-2) 

APEX 

18 

17.9 + 0.2 

3) 

*2CO(y = 3-2) 

JCMT 

14 

16.1+0.2 

4) 

*2CO(y = 4-3) 

APEX 

14 

16.4 + 0.5 

5) 

‘2C0(7 = 4-3) 

JCMT 

12 

6.7+ 0.2 

4) 

‘2C0(7 = 5-1) 

HIEF 

36 

1.2 + 0.02 

1) 

‘2CO(y = 9-8) 

HIEF 

20 

0.09 + 0.04 

1) 

*3co(y=i-o) 

IRAM 

22 

2.0 + 0.09 

1) 

*3co(y=i-o) 

SEST 

47 

0.45 + 0.02 

2) 

‘3co(y=2-i) 

SEST 

23 

1.3 + 0.04 

2) 

‘3C0(7 = 3-2) 

APEX 

19 

1.9 + 0.05 

6) 


Notes. The errors are Icr noise limits. 1) Danilovich et al. (20151, 
2)|01ofsson et al.| ( [T9^, 3)|De Beck et al.| (120101 ), 4) jSchoier et at] 
'2005| ), 5) this paper, 6) |Vlemmings et al. ( |2013| ) ^^^~The HIFI instrument 


on board the Herschel Space Observatory. 


Data for CO assembled from the literature are used to con¬ 
strain the circumstellar models. In the case of R Scl we have ten 
and four lines of ^^CO and ^^CO, respectively, in the radiative 
transfer analysis (Table |^, and in the case of S Set we have four 
*2C0 lines (Table [sjl. 

The CO line emission of detached-shell sources consists of at 
least two components, the emission from the present-day mass 
loss and the emission from the detached shell ( [Olofsson et al.j 
\\996) . These components are normally well separated in veloc¬ 
ity space since the velocity of the present-day wind is normally 
low 5 km s“^), i.e. about 30 % of the detached shell expansion 
velocity ( [Olofsson et al.[[1996[ ), e.g. in the case of S Set. How¬ 
ever, in the case of R Scl the situation is more complicated. The 
present-day wind velocity is 60% of the detached shell expan¬ 
sion velocity (a good estimate of the former is obtained from the 
^^CO(/ = 9 -8) line profile, ^ 8.5 km s“\ see Eig.[^. In addition, 
as shown by Maercker et al. (2012) using ALMA^CO(/ = 3-2) 
data, there is emission in the region between the present-day 
mass-loss envelope and the detached shell (in the form of a spi¬ 
ral pattern), and in the form of weak wings at velocities up to 
^19 km s“^ with respect to the systemic velocity. We have there¬ 
fore chosen to perform a fit only to the integrated (over velocity) 
line intensities in the ranges 0.75 Vq < \ P-Psys I < 1.05 Pg. The in¬ 
tegrated intensities in these ranges are listed in Table [^ The CO 
input data for S Set are given in Table 

In the Cl excitation analysis, we included the three ground- 
state levek_2Ro^_2^,_^^^ ^ 2 - The collisional data were taken 
from [Schroder et al.j ( [1991[ ) . 

The best-fit models were found using statistic defined as 


2 _ (4nod ^obs)^ 

^ 


( 2 ) 


where I is the integrated line intensities of the emission of the 
detached shell as observed (defined in Tables [^ and and in 
the models, N the number of observed lines, and cr the uncer¬ 
tainty of the intensity measurements, generally assumed to be 
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CTl 

O 



log(n[H2] / cm 

3.0 3.5 4.0 4.5 5.0 5.5 6.0 


D1 

O 



-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 


log(n[H2] / cm“^) 

3.0 3.5 4.0 4.5 5.0 5.5 6.0 


1 


8.4 

7.2 

6.0 


ai 

o 


2.4 




1.2 

0.0 



log(n[H2] / cm"^) 

3.0 3.5 4.0 4.5 5.0 5.5 6.0 




8.4 

7.2 

6.0 


4.8c§ 

3.6'^ 




2.4 

1.2 

0.0 


\og{M[H,]/M^) log(MlH2]/M0) log(MH2]/Mg) 

Fig. 2. (in logarithmic scale) for different abundances [10 ^ (left), 10 ^ (middle), 10 ^ (right)], where the density and kinetic 

temperature of the detached shell are free parameters, in the case of R Scl. 




V (km/s) V (km/s) 

Fig. 3. Observed (blue) and model (red) lines towards R Scl (A=APEX, H=HIFI, I=IRAM, J=JCMT, S=SEST). The left panel corresponds 
to a 3x10^ cm^/10 K model, and the right panel to a 2x10^ cm^/100 K model. In both cases fco = 10“^. The velocity scale is centred on the systemic 
velocity estimated to be -18.9 km s"k The dashed lines indicate the regions used in the fitting to the line intensities. 


Table 5. Integrated CO^line intensities for SSct in the velocity ranges results and vary the ‘^CO abundance, /nco, until a best fit is 
11 < I v-Vsys I < 19 km s obtained for a particular '^CO model. 


Line 

Tel. 

Beam 

["] 

[Kkrns-*] 

ref."^ 

'•^CO(7= 1-0) 

SEST 

45 

4.4 + 0.06 

1 

'2cO(7 = 2-1) 

SEST 

23 

3.0 + 0.07 

1 

'2CO(y = 3-2) 

APEX 

18 

0.56 + 0.1 

2 

'2CO(y = 3-2) 

SEST 

15 

1.0 + 0.09 

1 

Notes. The errors are Icr noise limits. 

1) Olofsson et al. ( 1996|), 


|De Beck et al.H2010] ) 


20 % unless otherwise noted. The reduced value is given by 
^red - I~ P)^ where p is the number of free parameters. 
In the determination of the detached-shell characteristics using 
the line data, the free parameters are the gas density and 
kinetic temperature. The abundance (with respect to H 2 ), 
/co, is also unknown, so we provide results for a number of 
models with different values for this parameter. In the Cl line 
emission modelling of the single observed line, the circumstellar 
model is fixed by the line modelling, and the only remain¬ 
ing free parameter is the Cl abundance with respect to H 2 , /ci- 
Similarly, for the ^^CO line modelling we use the ^^CO model 


5. Results 

5.1. Mass, density, and kinetic temperature of the detached 
shell 

In the case of R Scl, we have run models for five different CO 
abundances, fco = (0.01, 0.03, 0.1, 0.3, l)xl0“^. The last abun¬ 
dance applies to a carbon-rich medium where all the oxygen is 
locked up in CO (assuming a solar oxygen abundance of 5x10“"^ 
( [Grevess^ 2009| )), while the first applies to a medium where 
for example photodissociation has lead to a considerably low¬ 
ered CO abundance. The x^Qd ^^^ults are shown in Fig. 0 As 
can be clearly seen, reasonable fits (considering the uncertain¬ 
ties in extracting only the detached-shell emission) are found in 
large regions of the (7TH2,7'k)-plane. The solutions with the low- 
values are found along a curve from high-densities/low- 
temperatures to low-densities/high-temperatures (its position de¬ 
pends on the CO abundance), and acceptable solutions are found 
over broad ranges in density and temperature to the left of this 
curve. The high-density/low-temperature solutions correspond 
to optically thick line emission, and detached-shell masses that 
are much higher than acceptable (see Sect. |6.2[ ). For higher CO 
abundances, solutions of low-density/high-temperature models 
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H. Olofsson et al.: Cl line emission towards R Scl 


also provide good fits to the data. Results for some representative 
models are given in Table The predicted line profiles for one 
low-density/high-temperature model and one high-density/low- 
temperature model (both with fco = 10“^) are shown in Fig.|^ 

We have also determined the best-fit ^^CO/^^CO ratios for 
the presented models. The resulting ratios are given in Table 
and the predicted line profiles for one low-density/high- 
temperature model with fco = 10“^ is shown in Fig.|^ It is clear 
that the high-density/low-temperature models result in CO iso- 
topologue ratios that are unlikely. This is discussed further in 
Sect.lO 

One can conclude from this that with only the CO lines as 
constraints, and within the adopted homogeneous-shell model, it 
is difficult to estimate independently the density and kinetic tem¬ 
perature in the detached shell. The acceptable parameter space 
is large, but the low-density/high-temperature solutions, only 
present in the high CO abundance case, are the preferred option. 

We have performed the same analysis for S Set based on the 
^^CO input data given in Table We restrict ourselves to only 
two CO abundances in this case, 10“"^ and 10“^. The results are 
shown in Fig.[^ The general characteristics of the results are the 
same as those for R Scl, but the parameter space is more confined 
in this case, and the;^^^^ of the best-fit models are lower by about 
a factor of two. The latter is most likely an effect of the difficulty 
in isolating only the emission from the detached shell in the case 
of R Scl. 

5.2. CI/CO abundance ratio 

Based on the R Scl circumstellar model, constrained by the 
^^CO line data, we have estimated the abundance of neutral C 
compared to that of H 2 in its detached shell. It turns out that 
within a rather broad range of parameters (CO abundance, den¬ 
sity, and kinetic temperature) the CI/CO abundance ratio is rel¬ 
atively constant, and lies in the range 0.2-0.5 for the models 
with low to moderate optical depths in the ^^CO lines, i.e. the 
lower-density/higher-temperature models in TableJ^ The result¬ 
ing CI(^Fi - ^Pq) and the predicted CI(^F 2 - ^ 1 ) line pro¬ 
files are shown in Fig. |^for the 2x10^ cm“Vl00K//co = 10“^ 
model. The fit to the observed CI(^Fi - spectrum is good, 
indicating that most of the Cl line emission comes from the de¬ 
tached shell, but it is noteworthy that a not insignificant amount 
of CI(^Fi - ^Pq) emission may also come from the present-day 
slower wind. 

In the case of S Set we can only derive upper limits to the 
CI/CO abundance ratio, which are of the order of 1 (see Ta¬ 
ble 1^. This is not a particularly constraining limit, but it indi¬ 
cates that a substantial amount of the CO molecules in the de¬ 
tached shell are still not photodissociated unless the majority of 
the carbon atoms have been ionised. 


6. Discussion 

As has been shown, we can fit the observational CO data within a 
broad range of densities (and hence detached-shell masses) and 
temperatures. However, there are a number of arguments, pre¬ 
sented below in the case of R Scl, that will guide us towards the 
most reasonable description of the detached shell. 

6.1. Kinetic temperature of the detached shell 

The kinetic temperatures in “normal” circumstellar envelopes 
have been well constrained through the modelling of CO radio 



-20 -10 0 10 20 -20 -10 0 10 20 

V (km/s) 


Fig. 4. Observed (blue) and model (red) ^^CO lines towards R Scl 
(A=APEX, I=IRAM, S=SEST). The 2x10^ cmVlOOK and fco = IQ-^ 
model together with a ^^CO/^^CO ratio of 13 is used in calculating the 
model spectra. The velocity scale is centred on the systemic velocity 
estimated to be -18.9 km s"^. The dashed lines indicate the regions used 
in the fitting to the line intensities. 


lines ( |01ofsson]|2003| ). Using these results, we estimate the ki¬ 
netic temperatures at a distance of 10^^ cm from the star to be 
of the order 10 K. Thus, even though the temperatures in the de¬ 
tached shells may be elevated above this value, for example due 
to interaction with a previous wind, it appears unlikely that they 
are much warmer than about 100 K, which is also the preferred 
radiative transfer solution in the density/temperature space. 


6.2. Mass and line optical depths of the detached shell 


The detached-shell masses become very high for the higher- 
density/lower-temperature solutions, ^ 1 Mq. Although this is 
not impossible, it remains highly unlikely when considering 
for example the dust mass estimated for this detached shell, 
3x10“^ Me using dust-scat tered light (Gonzalez Delgado et al 


2003} Olofsson et al.||2010])^ and ^3 xl0~^ Me using dust ther- 
mal emission ( [Schoier et ^ 2005| . Adopting the 3x10“^ M© 
dust mass estiniate of [Schoier et al.| ( |2005| ) of the detached shell 
(based on thermal emission and therefore more reliable than 
those estimated from scattered light) and assuming a gas-to- 
dust mass ratio of 200 (applicable to the circumstellar medium), 
the detached-shell mass estimate is 6x10“^ Mq. Even though 
there are considerable uncertainties in this estimate, it points 
to the fact that a detached-shell mass above 0.1 M© is highly 
unlikely. In this limited range of the (nii^,T^)-pl 2 ine the lower- 
density/higher-temperature models, only present in the high CO 
abundance case, are preferred. 

The very high CO optical depths obtained for the high- 
density/low-temperature models are also unlikely since they re¬ 
quire ^^CO/^^CO abundance ratios in excess of 100, while the 
estima ted value for the ratio in the stellar photosphere 

is 19 ( [Lambert et al.|p^86| ). This value is not very far from 
the ^^CO/^'^CO line intensity ratios of 10 (see Table |^. Even 
though strong selective photodissociation of ^^CO cannot be 
excluded, the most reasonable interpretation is that the ^^CO 
optical depths are at most moderately high. Assuming LTE 
conditions, if we take as an indicative criterion that the opti¬ 
cal depth of the CO(/ = 3-2) is 1, this will happen at a radial 
H 2 column density of Nu 2 ~ 3x10^^ fco It should be 
noted that since the ^^CO results speak against high-density/low- 
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Table 6. Results for some representative models 


Source 

fco 

^^2 

Mh, 

T\dn 


'5CO/‘5CO 

CI/CO 



[cm“^] 

[Me] 

[K] 




RScl 

1x10-5 

2x10* 

3 

10 

3.2 

120 

0.056 


1x10-4 

7x105 

1 

10 

6.3 

390 

0.017 


1x10-5 

3x105 

0.04 

10 

3.9 

150 

0.045 


1x10-5 

2.4x105 

0.0035 

35 

3.7 

17 

0.22 


1x10-5 

1.5x105 

0.0022 

100 

2.4 

12 

0.35 

S Set 

1x10-4 

2x105 

0.07 

20 

1.6 


<0.7 


1x10-5 

2x105 

0.007 

100 

1.8 


<1.1 

2.6 

log(n[H 2 ] / cm"- 
2.8 3.0 3.2 

9 

3.4 3.6 



2.0 

log(n[H 2 ] / cm"9 

2.2 2.4 

2.6 




-2.0 -1.8 -1.6 -1.4 -1.2 
log(MlH2]/Me) 


-2.4 -2.2 -2.0 

logCMHJ/M^) 


I 




3.4 
3.0 
2.6 
2.2 
i.s| 

1.4 
1.0 
0.6 
0.2 


Fig. (in logarithmic scale) for different CO abundances [10 ^ (left), 10 ^ (right)], where the density and kinetic temperatures of the 

detached shell are free parameters, in the case of S Set. 
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Fig. 6. Observed (blue) and model (red) ^^CO lines towards S Set 
(A=APEX, S=SEST). The 2xl0^cmVl00K and fco = 10-^ model is 
used to calculate the model spectra. The velocity scale is centred on the 
systemic velocity estimated to be 14.9 km s"^. The dashed lines indicate 
the regions used in the fitting to the line intensities. The feature at the 
systemic velocity originates in the present-day wind and is not part of 
the modelling. 


Fig. 7. Left: The CI(^Pi - ^Pq) observed (blue) and modelled (red) 
line spectrum. Right: The CI(^P 2 “ modelled line spectrum. 
The 2xl02cmVl00K and /co = 10“^ 


model is used, together with 
fci/fco = 0-35, to calculate the model spectra. The velocity scale is cen¬ 
tred on the systemic velocity estimated to be -18.9 km s"^ The dashed 
lines indicate the regions used in the fitting to the line intensities. 


temperature models, they also speak in favour of a high CO 
abundance. 


6.3. Best-fit model for a homogenous detached shell 

Using the constraints on detached-shell mass and line optical 
depths, the models that provide the best fits to the line intensities 
are the lower-density/higher-temperature ones with a high CO 


abundance, 10“^. The best-fit values are = (2+l)xl0^ cm“^ 
(corresponding to Msheii = (2+l)xl0“^ M©) and Pk = lOO^yo^ K 
(with the caveat that a temperature above lOOK is unlikely). 
These result s are in excellent agreement with those of |Maercker| 
et al. ( 2015 ) based on ALMA ^^CO J= 1-0, 2-1, and 3-2 data. 
This corresponds to a mass-loss rate of 10“^ Moyr“^ during 
the ejection of the detached shell. The best-fit ^^CO/^^CO and 
CI/CO abundance ratios are 14+5 and 0.35+0.15, respectively. 
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H. Olofsson et al.: Cl line emission towards R Scl 


However, there are a number of other constraints provided 
by such processes as photodissociation and photoionisation, and 
we look at their consequences for our interpretation. 


6.4. Photodissociation of carbon species 

Our Cl results can, under certain assumptions, further constrain 
our modelling. Assuming that oxygen is of solar abundance and 
that all of it is initially locked up in CO, that all other carbon¬ 
bearing molecules are dissociated, and that there is no ionisation 
of carbon, we have for the ratio of neutral carbon and CO in the 
detached shell, 

M ^ 2[0/H]o(C/0-l) + /d/co,o 
fco (1 - /d)/cO,0 


where /d is the fraction of CO that is photodissociated, and /co,o 
is the abundance of CO assuming all oxygen is locked up in CO 
(i.e. 10“^). Adopting the value C/0 = 1.34 for R Scl from 


Lam- 


bertetal.| ( [T98^ we find fci/fco = 0.34 if /d = 0 (i.e. fco = fco.o). 


which is in excellent agreement with our estimated ratio for the 
best-fit model. This may be fortuitous, but may, on the other 
hand, indicate that our assumptions are basically correct. Our 
results are consistent with the predictions as long as /d^O.l, 
i.e. only a small fraction of the CO molecules are photodisso¬ 
ciated. On the other hand, if carbon ionisation is effective, we 
get no co nstraint on the CO abundance in the detached shell (see 
Sect |6.5| ) and the fact that the CI/CO abundance ratio is 0.35 is 
fortuitous. 

We can make an estimate of the effectiveness of CO 
photodissociati on. Its un-shielde d photodissociation rate 


IS 


2.6x10“^^ s“^ ( Visser et~^|2009 l, but - since it is dissociated 
in lines - it is very efficiently self-shielding and in addition H 2 
lines contribute to the shielding. [Visser et~SI] ( |2009| ) have cal¬ 
culated the extra shielding provided by CO itself and H 2 , and 
as a guideline CO and H 2 column densities of 10^^ cm“^ and 


10^^ cm“"^ (applicable to our best-fit model), respectively, gives 
an extra shielding by a factor of 25 (although the parameters 
used in the calculation of this do not match those of our detached 
shell), meaning a life time of the order of 3000 years, i.e. longer 
than the age of the R Scl detached shell. This makes it likely that 
the CO molecules have survived in the detached shell. 

With the assumptions outlined we have shown that the CO in 
the detached shell of R Scl is affected very little by photodisso¬ 
ciation (which is also consistent with the fact that detached CO 
shells of ages approaching 10"^ years, e.g. TT Cyg and S Set, are 
found, i.e. the CO molecules are still present at these ages). The 
carbon atoms must therefore come from the full dissociation of 
the main carbon-bearing species, apart from CO, and whatever 
carbon that remained in atomic form after the formation of CO 
and other molecules. Of course, the propensity for photodissoci¬ 
ation is strongly dependent on the dumpiness of the circumstel- 
lar medium and this may very well differ significantly from one 
detached shell to another. 


6.5. Photoionisation ofC 

The strength of the observed Cl line is consistent with a high CO 
abundance. However, this conclusion requires that the amount of 
carbon that has been ionised is very limited. 

The observational constraint on the amount of ionised carbon 
is restricted to a spectrum of the CII line at 158 pm in the ISO 
archive, which is difficult to assess owing to the poor and limited 
spectral resolution. It corresponds to an integrated line flux of 


10“^^ W cm“^, which we will treat as an upper limit. Because of 
the high energy of the upper level (91 K), we can expect to see 
CII line emission only in the higher-temperature models. Using 
our best-fit model, and assuming conservatively that all the CII 
line emission comes from the detached shell, this results in an 
upper limit to the CII/CI abundance ratio of 7 in the detached 
shell, i.e. not a conclusive upper limit in this context. 

We can also make an estimate of the efficiency of carbon 
photoionisation. The ionisation rate of carbon is esti mated to 
be 3x10“^^ s“^ (for a normal interstellar radiation field: {Draine 


|1978j )] and (van Dishoeckj 19881 )). With no extinction, the ionisa¬ 
tion time scale is only 150 yr. This means that any un-shielded 
carbon should have been ionised long before it reached the loca¬ 
tion of the detached shell. We can make an estimate of the shield¬ 
ing effect. A radial H 2 column density of 10^^ cm“^ through the 
detached shell corresponds to an extinction Ay 


1 (Bohlin et^ 

|1978| ) (assuming interstella r conditions). The d epth correction 
is therefore about ^ 30 ( van Dishoeck 1988| ), increasing the 
time scale to 3000 years, i.e. longer than the age of the R Scl 
detached shell. However, the best-fit model for a homogenous 
detached shell has a radial H 2 column density of only 10^^ cm"^, 
i.e. not enough to effectively protect the carbon against ionisa¬ 
tion. 


6.6. The ^^COA^CO ratio in the detached shell and its 
implications 

There is one more constraint to be discussed, the CO isotopo- 
logue ratio. In their study of the ^^CO/^^CO ratio, based on 
ALMA images which resolve the emissions, in the circumstel- 
lar environment of R Scl, jVlemmings et al.j (2013| ) found a 
^^CO/^^CO ratio of 60 close to the star and (on average) 20 
in the detached shell. They interpreted this in terms of selec¬ 
tive dissociation of close to the star and a re-formation 
of ^^CO, through the chemical fractionation reaction, in the de¬ 
tached shell. This requires that the detached-shell temperature is 
low, especially lower than the difference in ground-state energy 
between the two isotopologues, which in the case of and 
is 35 K. Thus, these data, if correctly interpreted, advocate 
a lower-temperature solution. 


6 . 7 . Clumpy medium 

What we have discussed so far cannot be fully reconciled if we 
assume a homogenous detached shell. In particular, protection 
against photoionisation of C, and to some extent also the pho¬ 
todissociation of CO, require higher H 2 column densities. This 
may point towards a clumpy medium, where the connection be¬ 
tween density, shell width, and mass is offset, and for example 
higher-density/lower-temperature solutions for a high CO abun¬ 
dance and a limited detached-shell mass become possible. A 
clumpy medium is also consistent with the small-scale morphol¬ 
ogy of the detached shells observed in sufficient detail (Olofssonj 
jet al.|2000l|MT0l|Maercker et al.|2012D . 


Bergman et al.j (1993| ) and jOlofsson et al.j (1996| ) have pre¬ 


sented a model of a clumpy detached shell. The radiative trans¬ 
fer of CO, HCN, CN, Cl and CII is solved simultaneously with 
the physical evolution of a clump as it moves away from the 
star. The photodissociation of the molecular species and the 
photoi onisation of C are taken into account. jOlofsson et al. 
( 1996 ) present the results for a 10“^ M© clump of original radius 
10^^ cm and temperature 750 K. The evolution is followed over 
1.5x10'^ yr. In 2300 years the clump has expanded to a radius of 
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4x10^^ cm, has an H 2 density of J^2xl0^cm“^, an H 2 column 
density of 1.5x10^^ cm“^, and a temperature of 35 K. It was 
further shown that a thousand such clumps in a detached shell 
gives a very good fit to both the and line intensity 
data towards R Scl using an abundance ratio ^^CO/^^CO = 19 
(with a much broader shell than the present observational data 
indicate, but it is expected that this has only marginal effects on 
the line intensities). 

Initially, the Cl emission comes from carbon released from 
HCN (via CN), in this specific case at an abundance 125 times 
lower than that of CO. Only after 5000 years is there a con¬ 
tribution from photodissociation of CO, i.e. well beyond the age 
of the R Scl detached shell. The predicted line intensity ratio 
between the Cl line and the CO(/ = 3-2) line at 2300 years is 
0.05 (scaling their result to the Cl abundance derived by us and 
using the APEX beams at the relevant frequencies), remarkably 
close to the observed line intensity ratio of 0.03 at APEX. Ei- 
nally, the estimated CII to Cl abundance ratio is 0.5, i.e. 30 % 
of the carbon has been photoionised. 

This makes it likely that a clumpy medium will provide an 
explanation for the observational data, even though we have not 
made any fine-tuned modelling of such a medium that takes into 
account for example the much smaller width of the detached 
shell as shown by the more recent observations. 


7. Conclusions 

The physical properties of detached CO shells around carbon 
stars are difficult to constrain even when using a substantial num¬ 
ber of CO (and ^^CO) lines obtained with single-dish telescopes. 
It has also proven very difficult to detect any other molecular 
line emission from these detached shells. Our detection of the 
CI(^Pi - ^Pq) line is a step forward and provides additional infor¬ 
mation. The best fit to the observed CO line intensities, assum¬ 
ing a homogeneous detached shell, is obtained for a detached- 
shell mass of ^ 0.002 Mq and a temperature of 100 K if the 
CO abundance with respect to H 2 is 10“^ (the ^^CO line emission 
excludes models with low CO abundances). The CI(^Pi - ^Pq) 
line is detected at an intensity that results in a CI/CO abundance 
ratio of 0.3 in the best-fit model. This is consistent with the car¬ 
bon coming from the full photodissociation of all carbon-bearing 
species except CO, i.e. a high CO abundance is also inferred in 
this way assuming that the ionisation of carbon is very limited. 
An ISO spectrum of the CII 158 pm line gives an upper limit to 
the CII/CI abundance ratio of 7 in the case of our best-fit model, 
i.e. not a conclusive upper limit in this context. It would therefore 
be worthwhile to redo these observations of the CII line at much 
higher sensitivity. It is further concluded that a situation that ful¬ 
fils all the observational constraints, e.g. the required column 
densities to limit CO photodissociation and C photoionisation, 
only exists in a clumpy medium. However, the details of such 
a medium (e.g. the size distribution, density, and temperature of 
clumps) remain uncertain. Substantial progress can only be made 
through detailed observations of the small-scale structure of the 
circumstellar medium using instruments such as ALMA, where 
the properties of individual clumps are determined. 
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